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Abstract 

Recently it has been discovered that some nonlinear evolution equations in 2 + 1 
dimensions, which are integrable by the use of the Spectral Transform, admit lo- 
calized (in the space) soliton solutions. This article briefly reviews some of the 
main results obtained in the last five years thanks to the renewed interest in soliton 
theory due to this discovery. The theoretical tools needed to understand the un- 
expected richness of behaviour of multidimensional localized solitons during their 
mutual scattering are furnished. Analogies and especially discrepancies with the 
unidimensional case are stressed. 



1 Introduction 

1.1 Special features of solitons in two dimensions 

Since the discovery of the soliton in 1965 by Zabusky and Kruskal a large new domain of 
mathematical physics developed and is believed to have reached maturity. It is generically 
called the soliton theory. Its principal mathematical tool is the so called Spectral (or 
Scattering) Transform that is used to solve a large class of nonlinear evolution equations 
in 1-1-1 (one spatial and one temporal) dimensions. 

Some of these equations, in particular the nonlinear Schrodinger (NLS) equation and 
its generalizations, can be obtained in an appropriate multiscale limit from a very large 
class of nonlinear dispersive equations. Therefore, it is not astonishing that applications 
of soliton theory are percolating through the whole of physics, especially quantum field 
theory, solid state physics, nonlinear optics, plasma physics and hydrodynamics, and 
other natural sciences. 

The most impressive phenomenon in the theory and in the applications is the existence 
of solitons, i.e. (localized) coherent structures that mutually interact preserving their 
individuality. 

In the last decade many efforts have been made to extend the soliton theory to non- 
linear evolution equations in 2-1-1 (two spatial and one temporal) dimensions. In fact 
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the Spectral Transform was extended to dispersive nonlinear evolution equations in 2+1 
dimensions but it was generally admitted the lack of two dimensional localized solitons. 
Only recently, in 1988, it has been discovered (by Boiti, Leon, Martina and Pempinelli) 
that all the equations in the hierarchy related to the Zakharov-Shabat hyperbolic spectral 
problem in the plane have (exponentially) localized soliton solutions. The most repre- 
sentative equation in the hierarchy is the Davey-Stewartson I equation (DSI), which 
provides a two dimensional generalization of the NLS equation. 

This discovery has stimulated a renewed interest in soliton theory. The first results 
are very promising. In particular soliton solutions display a richer phenomenology than in 
1+1 dimensions. This opens the way to applications in multidimensions, which, hopefully, 
are expected to be even more interesting than in one space dimension. 

In contrast with the 1+1 dimensional case the time evolution of the solution of the 
Davey-Stewartson equation is not uniquely determined by the initial data. In addition 
one has to give, at all times, boundary data. If they are chosen to be identically zero 
solitons cannot be present. For a convenient choice the solution can contain solitons but 
not necessarily does. In the affirmative case the boundary data fix the kinematics of the 
incoming and outgoing scattering solitons, i.e. their velocities and locations in the plane 
in the far past and in the far future. The initial data fix the dynamics of the interaction. 

The scattering of the solitons can be inelastic and they can change shape and also 
exchange mass (energy or charge according to the specific physical interpretation). In 
fact, while the total mass of solitons is conserved, the mass of the single soliton, in 
general, is not preserved by the interaction and solitons can also simulate inelastic scat- 
tering processes of quantum particles as creation and annihilation, fusion and fission, and 
interaction with virtual particles. 

1.2 Guidelines for additional reading 

One main feature of the nonlinear dynamical systems is that different non equivalent 
approaches are possible and each one is useful and clarifying from a special point of 
view. 

Also multidimensional solitons can and, effectively, have been studied by using dif- 
ferent tools. In this article, mainly for lack of space, we made the (questionable) choice 
of using only the Backlund transformations and a special version of the Spectral Trans- 
form. In fact, we collected, reorganized and simplified the results that the reader can find 
scattered, with some minor additional details, in references In the first article of 

the list the localized solitons in the plane have been discovered and in the second one a 
preliminary analysis of the properties of the new Spectral Transform needed to describe 
them is performed. These two papers opened the way to a deeper understanding of inte- 
grable nonlinear evolution equations in multidimensions and to the search and discovery 
of other dynamical systems admitting localized coherent structures. 

Another special version of the Spectral Transform theory different from that one 
presented in this paper has been developed in ||]-||ll|. This alternative theory is relevant 
because it has been used to clarify the role played by the boundaries and to show that 
multidimensional solitons, in contrast with the one dimensional solitons, can interact 
inelastically. These authors suggested to call these solitons dromions in order to stress 
the fact that they can be driven everywhere in the plane along tracks (dromos in greek) 
by choosing a suitable motion of the boundaries. 
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In both versions the Spectral Transform for the nonstationary Schrodinger equation 
(with a potential vanishing at large distances in every direction except a finite num- 
ber) plays a fundamental role. To extend the theory, originally developed for potentials 
vanishing at large distances p^ , to this case it has been necessary to introduce a new 
mathematical entity resembling the resolvent of the linear operator theory. For lack 
of space this topic has been skipped in this article. The interested reader can consult 
Ip^ , [l^ . The multi-soliton solution of the Kadomtsev-Petviashvili I (KPI) equation, 
that is needed in order to build the multi-soliton solution of the DSI and DSIII systems, 
has been derived in and in its most general form in . 

Also other methods have been used to get the multidimensional solitons, the bilinear 
approach , the quantum machinery of creation and annihilation operators |9| , the 
Grammian determinants [ po[ , the Darboux transformations and the dressing method 
|p2| . The most general form of the multi-soliton solution together with the remark that 
the number of solitons is not necessarily conserved were first given in . 

The search of multidimensional solitons has been extended by different authors to 
other nonlinear evolution equations in 2 + 1 dimensions p^-p6|. Some interesting at- 
tempts using the 9-method and a direct method have been done also in higher dimen- 
sions, precisely in |2^] and 

In |Q-|3^ it has been shown that the DSI and, successively in that also the 
DSIII system can be obtained by using a multiscale limit starting from a very large 
class of dispersive nonlinear equations. In particular the DSI equation with boundaries 
compatible with solitons can be embedded into the KPI equation |3^. These results 
open the way to the search of physical applications other than hydrodynamics, studied 
in p8| , |39| , and references quoted therein] . In this respect a more discouraging analysis 
is made in [Q. This last paper, mainly dedicated to the consequences of the lack of 
conserved quantities for DSI, contains also an interesting examination of the literature 
dealing with the discovery of the DSI system and with the early attempts to build the 
corresponding Spectral Transform. 

For the Hamiltonian version of the DSI system and its quantum extension see [^-jisll. 
A Hamiltonian version of the DSIII system appears already in and a bi-Hamiltonian 
version in |^7| . For the general method for building intcgrable nonlinear evolution equa- 
tions starting from a 5-problem that allowed to prove that DSIII is intcgrable see psj - 
pO | and, specifically, For more details on the localized soliton solutions of DSIII see 
|p2 |. For the DSII system which is related to the elliptic version of the Zakharov-Shabat 
problem in the plane see Q for the Spectral Transform theory and |5^ for the 
singular soliton solutions. 

Finally, between the different excellent existing books on solitons let us suggest to the 
reader p7| and , as the most accurate and complete for the nonlinear evolution equa- 
tions associated to the Schrodinger and to the Zakharov-Shabat spectral equation in 1+1 
dimensions, respectively, and ||59t| as the most updated and comprehensive. In particular 
for those who want to have a general overlooking on the subject and a rich bibliography 
to pick over this book is particularly recommended. A new book on multidimensional 
soliton has been just announced jGOt . For the Backlund and Darboux transformations 
see 1^ and |^ . Special attention to the algebraic and geometric approach to the soliton 
equations has been dedicated in and to the multidimensional case in [ |6^ . 

In all these books the reader can find many useful references for the problems consid- 
ered in this paper and for related problems. Here we want to quote only some references 
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more relevant historically or more close to our specific approach. Precisely [^5[ |6^ for 
the discovery of the one dimensional solitons, ^ for the extension of the theory 
to the nonlinear Schrodinger equation, |69 for the introduction of the so-called d- 



method, 0-j7|| J5j, |]J, [[75|~|82l for the extension of the Spectral Transform to 2 + 1 
dimensions, |76|]-||79[ in particular for the introduction of the "weak" Lax representation 



of the integrable equations in 2 + 1 dimensions, [^l, 82 for the Backlund transformations 
and for the extension of the Backlund transformations to 2 + 1 dimensions. 

By offering these guidelines for recovering interesting references we hope to redress 
the omissions due to oversight that have certainly occurred. Anyway, we apologize for 
this to the reader as well to the authors who may have unjustifiably excluded. 

1.3 The Davey-Stewartson I equation 

The Davey-Stewartson (DS) systems model the evolution of weakly nonlinear water 
waves that travel predominantly in one direction, are nearly monochromatic and are 
slowly modulated in the two horizontal directions. We are interested in the special DS 
system (DSI equation) that one gets in the shallow water limit when the effects of the 
surface tension are important. In characteristic coordinates and dimensionless form the 
DSI equation is a system of two coupled equations 

iqt + quu + qvv - iVu + Vv- cra\qf)q = (1.1) 
2(puv ^ cro{\q\'^)u + (7o{\q\'^)y, CTq = ±1 

where q{u,v,t) is the (complex) envelope of the free surface of the water wave we are 
considering and (/?(«, v, t) is the (real) velocity potential of the mean motion generated 
by the surface wave. 

It is worth to stress that the DSI system is not necessarily placed in the context of 
water waves. Indeed, it has been shown that a very large class of nonlinear dispersive 
equations in 2+1 dimensions reduces in an appropriate asymptotic limit to the DSI 
equation and therefore we expect it to arise in many different physical applications. 

To exhibit explicitly the bound ary value of f at large distance in the {u, v) plane 



allowed by the second equation in (LI) it is convenient to introduce the two fields 



= -(^^ + iao|# (1.2) 
A^'^=^u-\<r,\q\' 
and to rewrite the DSI equation as 

*gt + gun + g™ + (A(i) -^(2))^^0 (1.3) 



1 



f dv'a^{\q\^)u + a^i^\u,t) 
where a(^^ and Oq^^) are the arbitrary boundaries. 
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In fact we could make another not equivalent choice and write the DSl equation as 
iqt + quu + qvv + {A^^^ - ^(2))^ = (1.4) 

A(^^=-- (r + r \du'ao{\qW + A^^\v,t) 
A^^^^Ur + r )dv'ao{\qnu+A^o\^>t) 

^ \j —oo J-\-oo/ 

where now Al^\v,t) and A^^\u,t) are the arbitrary boundaries. 

It can be shown that the proper boundary conditions to be chosen are dictated by 
the specific multiscale limit one is choosing in getting the DSI equation. Specifically, 



equation (1.3) can be obtained via a multiscale limit from the Kadomtsev-Petviashvili 
(KPl) evolution equation while maintaining well posedness in time. 

In order to solve the DSI equation it is convenient to introduce its more general two 
component version 

) + [A,Q]^Q (1.5) 



where 



and 



with 



or 



= r du'iQ\+a^^\v,t) (1.8) 



AW=--(/ +/ )du'{Q'),+4'\v,t) (1.9) 
4 



+ / ) dv' {QX + Al^\u,t) 
J +00/ 



according to the boundary conditions we choose. The previous considered equations are 
simply obtained for 

r = aoq (1.10) 
(where q denotes the complex conjugate of q) and are called reduced DSI equations in 



contrast with (1.5), (1.8) and ( |l.9| ) which can be named DSI equations with no additional 
specification. 

The DSI equation can be obtained as the compatibility condition (Lax representation) 

[Ti,r2] = o (1.11) 
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for two underlying linear operators 



TiiQ) ^ + <73dy + Q 

T2[Q) = idt + azdl + Qdy - iasQ:. + \Qy + A 



where 



1/ 

2 + 



y 



\{u-v). 



(1.12) 
(1.13) 

(1.14) 



and (73 is the Pauli matrix. 

The first of these operators is the Zakharov-Shabat hyperbolic spectral operator in 
the plane and can be considered to define a linear spectral problem in which Q plays the 
role of the data. In the case ( |l.8| ) with o'q^ and Og^^ considered as boundaries given at 
all times one can define the Spectral Transform of Q and solve the initial value problem 
for the DSL This is explicitly done in section 3. 

In the case (|l.9|) one can introduce the Hamiltonian 



H 



dudv 



r(dl + dl)q ~ \qv(dudZ^ + d,dZ^)qr + {Af - J^^^)qT 



and the canonical Poisson brackets 

'^F8G SFSG] 
Sq 6r Sr 6q 

where q and r are the conjugate variables. Then the equations of motion 

{q,H}, rt = {r,H} 



{F,G} =ijjdudv 



qt 



(1.15) 



(1.16) 



(1.17) 



yield the DSI equation. The problem of defining a Spectral Transform is completely 



open and, moreover, only in the special case A, 



(1) 



A 



(2) 



it has been shown that 



the DSI equation is Hamiltonian with a continuous infinity of independent commuting 
constants and is completely integrable in the Hamiltonian sense. This case is usually 
named Hamiltonian case. 



1.4 The Davey— Stewartson III equation 

There is another nonlinear evolution equation that can be associated to the Zakharov- 
Shabat hyperbolic spectral operator in the plane and that admits localized soliton solu- 
tions. 

To get it we need to introduce a weaker form of the Lax representation (1.11). Pre- 
cisely, we search for a second spectral operator T2 that commute with the Zakharov- 
Shabat spectral operator Ti in (1.12) only on the subspace of the eigenfunctions of Ti 
("weak" Lax representation) 

Ti7/. = 0, [ri,T2]V' = 0. (1.18) 

For 

Ti = 2diag(9„,9,) + Q (1.19) 
T,^^^^ + ^l + ^',+A+f^ '^^ ) (1.20) 
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we have the so-called DSIII system 



,) + [A,Q]^0 (1.21) 

(1.22) 

(1.23) 

to be compared with the DSI system i n (|l.5D and (1.8) or ( |1.9| ). The DSIII system, as 
DSI, is compatible with the reduction (1.10). 

Of course also DSI can be obtained by using a "weak" Lax representation instead of 
the usual "strong" one in ( pTlll ). Then, one can choose for T2 instead of ( |1.13| ) 



iQt - 


- 0-3 {Qvv — 











T2 = idt + dl - dl 



A 



-qu 




(1.24) 



and get again DSI. The theory of the Backlund transformation and of the Spectral 
Transform developed in the following sections does not ch ange if the st artin g Lax pair is 
"weak" or "strong". Therefore, comparison of definition ( |l.2(]| ) with ( 1.24 ) makes clear 
that one cannot expect any difference, apart some signs, between formulae for DSI and 
DSIII. 

In particular the time evolution of i/' for DSIII has to be fixed as follows 



T2V = -k^^ 



to be compared with ( 3.62 ) for DSI. 

The DSIII system admits for convenient boundaries of the form 



Ad) 



du' (g2 



,(1) 



{v,t) 



(1.25) 



(1.26) 



dv' {Q'')u + ai\u,t) 



localized solitons of the same shape of those of DSI, which exhibit similar dynamical 
phenomena but evolve differently in time. 



We write here the one soliton solution 
2Acj77 exp [iip] 



D 



2/iQpexp[— i(/3] 
D 



where 



LP = 

D 
6 



^sRU + AsRW + (A| - A4 - A^i 
= 27(cosh^i + cosh^2) + exp(^2), 7 
-- -/igii - Aau + 2(A3A5R - HQHiit)t, 
-- /Ltcju - Aqu + 2(AcjA!ff + AiaMa?)* 



(1.27) 



(1.28) 
(1.29) 
(1.30) 
(1.31) 



to be compared with ( ^.40 )-(2.44). As in the DSI case the complex parameters A = 
AsR + iX^ and /i = /ij} + ifi^ are the discrete eigenvalues of the associated Zakharov- 
Shabat spectral problem and p and 77 are arbitrary complex constants satisfying the 
conditions 7 G H and 7(1 + 7) > 0. The multi-soliton solution can be easily obtained 
by taking the same lines we choose in the following sections for the DSI equation. 
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2 Solitons via Backlund Transformations 



The Backlund transformations have their origin in work by Backlund in the late nine- 
teenth century and are, therefore, the oldest tool used in exploring nonlinear integrable 
systems. Much more recently many different sophisticated and powerful methods have 
been developed, in particular the Spectral Transform and the dressing method. How- 
ever, in our opinion, the Backlund transformations remain the simplest way for getting 
the soliton solutions. Moreover, because, under appropriate circumstances, a reiterated 
application of the Backlund transformations generate a sequence of solutions by a purely 
algebraic superposition principle they can be used to study the interaction properties of 
the solitons. 

The simplest way to generate a Backlund transformation is to use the gauge invari- 
ance of the linear spectral problem associated to the nonlinear evolution equation one is 
considering. The gauge that generates the Backlund transformation is called Backlund 
gauge. The localized soliton solutions of the DSI equation, with boundaries of the form 
in (l.S), were for the first time derived by using these special gauge transformations. 
Successively they have been rederived by using the techniques of the Spectral Transform. 
But, in the case of the Hamiltonian DSI equation, we have not, presently, at our dis- 
posal the Spectral Transform or the dressing method and, consequently, in order to get 
explicit solutions we are left with the necessity to generalize the Backlund gauges such 
as to include also the special form of the boundaries in ( |1.9| ). To choosing these bound- 
aries corresponds to imposing to the solutions nonlinear constraints, which can be solved 
only by using the additional freedom at our disposal in the generalized Backlund gauges. 
We are able to write explicitly infinite wave solutions with constant and periodically 
modulated amplitudes. 



2.1 Generalized Backlund gauge transformation 

Once given a solution Q of the DSI equation we want to generate a new solution Q' of 
the same equation by using a convenient gauge operator B that transforms according to 
the equation 

i:' = B{Q',Q)^ (2.1) 
the matrix solution ip of the principal spectral equation 

Ti(Q)V' = (2.2) 

for Q to the matrix solution ip' of the same spectral problem 

Ti(g')V'' = (2.3) 

for Q'. 

It is easy to verify that if B satisfies 

Ti{Q')B{Q', Q) - B{Q', Q)T^{Q) - (2.4) 

T2{Q')B{Q', Q) - B{Q', Q)T2{Q) = (2.5) 

then Ti{Q') and 12 (Q') satisfy the same compatibility condition 

[ri,r2] = o (2.6) 
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as Ti{Q) and T2{Q) and, consequently, Q' satisfies the DSI equation. In this case the 
gauge B is cahed Backlund gauge and the equations (2.4) and (2.5) yield, respectively, 
the so-called space and time component of the Backlund transformation. 

Non trivial Backlund gauges are polynomial in the operator dy. We are interested in 
the most general Backlund gauge of first order of the form 



B{Q',Q)^ady+Bo{Q',Q) 



(2.7) 



with a a constant diagonal matrix and Bq a matrix. By inserting it in (2.4) we get its 
functional form 



B(g', Q) = ady - \a^{Q'a ~ aQ) ~ \(JiaI{Q'^ - 

and the space component of the Backlund transformation 

Q' [/3 - iaag X(Q'2 - Q2)] _ _ 1 ^,^3 X(Q'2 - Q'^)\ Q 
~ \a^{Q'a - aQ)^ - \{Q' a + aQ)y = 0. 

The matrix operator X is defined by 

T={d,,+a^dy)-^ 
and the diagonal matrix /3 is subjected to the constraint 

(d^ + G^dy)l5 = Q, 



i.e. it is of the form 



where 



/3i(v,i) 
h{u,t) 



V 



(2.8) 

(2.9) 

(2.10) 

(2.11) 
(2.12) 
(2.13) 



and /3i and (5i are arbitrary functions. Note that, in contrast with the 1+1 dimensional 
case, the matrix /3 that plays the role of 'constant of integration' in the solution of (2.4) 
admits also a space dependence. This additional freedom will be used in the following 
for getting soliton solutions of the Hamiltonian DSI equation. 

By inserting B{Q', Q) in ( pj] ) we get the time component of the Backlund transfor- 
mation, which can be shown by use of (2.9) to be equivalent to the DSI for Q' , and two 
additional partial differential equations 



[P ~ iaa3X(Q'2 - 0^)] + \aa^{^ - A) + \a{Q'^ - g') = 



(2.14) 



\ {a - \aa^ X(g'2 - g2)] ^ + (A' - A) [/3 - \aG^ I{Q'^ - Q^)] 
-\a{A' + A)y + \acy^{Q'^ + Q^)^ - \ [(9, + a^dy)Q'\ Q'a 
-\aQ[{d^ ~ a^dy)Q\ + \{d^ - (T^dy){Q' aQ) = 0. 



(2.15) 



These two equations can be used for determining the field A' and the admissible /3's. It 
can be verified that they are compatible with the equation (2.11) for (3 by applying the 
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operator (dx + crsdy) to both of them and by showing that the two obtained equations 
are identically satisfied for Q and Q' solutions of the DSI equation and of the space 
component of the Backlund transformation. 

The Backlund gauge B{Q', Q) for general a and (3 in (^) can be obtained by com- 
posing two simpler Backlund gauges that are called elementary Backlund gauges of the 
first and second kind. They are obtained by choosing, respectively, 

and 

"-=(2 ^-=(i MM))' ^'-''^ 

and are noted Bj{Q' , Q; A) and Bji{Q' , Q: ii). They are not compatible with the reduc- 
tion 

r^aoq (2.18) 

but by composing two elementary Backlund gauges of different kind one can get finally 
a solution satisfying the reduction. 

This procedure simplifies radically the computation because, according to the gen- 
eral feature of the Backlund transformations, the recursive application of the Backlund 
transformations can be achieved by purely algebraic means. An additional simplification 
is obtained by imposing the commutativity of the diagram 

Qi 

Q' 

Biiitj\ / Bi{X) 
Qii 

which represents symbolically the following equation 

B{Q', Q- A, li) = BiiiQ', Qi- tJi)Bi{Qi, Q: A) 

^ Bj{Q',Qii;X)Bji{Qji,Q;fi) (2.19) 

where B{Q' , Q; A, /i) is the Backlund gauge in ( ^.^ ) with 
2.2 Localized soliton solutions 



The most natural choice for the parameters entering ( 2.16| ) and ( 2.17 ) is of course to take 
A and /i constants, say 

\{v,t)=i\ fi{u,t)^ifi, (2.21) 



where A and ^ are complex constants. Then we choose the starting solution Q of (1.5) 
to be zero together with its related auxiliary field A (note, however, that the vanishing 
of Q does not imply the vanishing of A). 
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With these choices the commutation relation (2.19) becomes 

B(Q, 0; A, ^i) = Bjj{Q, Q,; /i)B,(Qj, 0; A) = Bj{Q, Qjj: X)Bji{Qjj,0; /i). (2.22) 
By equating to zero the coefficients of the powers of dy it results that 



where 



rj = p{u, t) exp(— iAw), qu = r]{v, t) exp{i^u), (2.24) 

q=—— , r = (2.25) 

1 + jTiqii 1 + jriqii 

and that the operator I satisfies the equation 

AQ^)^\q{Qi~Qii)- (2.26) 

Here above p and rj are arbitrary functions of the space variable, that can be represented 
by the following transforms 

p{u,t)= 1 1 dkAdke-'''''p{k,t), (2.27) 



r]{v, t)^ J J dkAdk e'''''fi{k, t), (2.28) 
with the explicit time dependencies 

p{k, t) = p{k, 0) exp[i(fc2 + A2)t], (2.29) 

fl{k, t) = fi{k, 0) exp[-i(fc2 + p^)t], (2.30) 

where /5(fc, 0) and fi{k,0) are arbitrary. 

Everything now lies in the choice of these two arbitrary functions. We prove hereafter 

that the requirement that the solution Q is localized and obeys the reduction condition 



(1.10) determines uniquely the structure of the arbitrary functions p and fj. To that end 
we redefine 

p(u,0) = 2e"'^"S'(a), (2.31) 

?7(i;,0) = 2e^-^"r(T), (2.32) 

with the new variables 

— ^ 

(7 = r exp[i(/i — /i)u], (2.33) 

p- p 

- exp[-i(A- A)t;]. (2.34) 



A- A 



In these variables, the reduction condition (1.10) takes the following form 



[1 + 5(a)f (r)]^ = [1 + 5(a)T(r)]^, (2.35) 
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which can be solved exphcitly and admits 4 independent solutions. The requirement of 
the localization selects the only solution 



S = aa + b, 
T — ar + c, 



ha — ba — 0, 
ca — ca — 0, 



(2.36) 



where a, 6 and c are complex constants. Then, by means of and we arrive 
at 

P{k,Q) ^ p[5{k - ^l) + 5{k - (2.37) 

77(fc,0) = 7/[(5(fc- A) + (5(fc-A)]. (2.38) 

Here above "q and p are arbitrary constants. The requirement of the reduction condi- 
tion gives the following constraint on these constants 

p{p-ii)^aon{^~l)- (2.39) 
Therefore the localized one-soliton solution to the system (|l.5| ) (1.8) can be written 
2Aa77exp[i(^] 2/icjpcxp[— 



D ' D 

with the following definitions 

if = /isRM + AsRW + (A| - A| + /4 - Ai^)*, 
D = 27(cosh^i + cosh^2) + exp(^2), 7 = 

= - -^3" + 2(A3A5ff + p,c3p,u)t, 

6 = - + 2(A3A3? - /^S/iSR)i, 

and where Asr denotes the real part of A and Aq its imaginary part. 

For Aa/iQ ^ and 7(1 + 7) > 0, the above formulae describe a two-dimensional bell- 
shaped envelope of the plane wave exp(±z</j), exponentially decreasing in all directions 
in the (m, u)-plane and moving without deformation at the constant velocity 



(2.40) 

(2.41) 

(2.42) 
(2.43) 
(2.44) 



(2.45) 



The initial position of the soliton is arbitrary, in other words we may translate the 
space variables according to the transformation 



u - uo, 



V - Vo, 



uo,vo e M, 



(2.46) 



and, consequently, the soliton is defined by means of eight real parameters. 

By inserting the value of T {Q^) obtained in ( 2.26 ) into ( 2.14 ), rewritten for the gauge 
B{Q, O; A, /x) derived in ( 2.22| ), we get the auxiliary field A 



(2.47) 
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Finally, from ( pTl| ) by using the gauge operator B{Q,0; X, fi) we obtain the following 
eigenfunction relative to the one soliton solution 

/ rig 2q \ 

^{x,y,k)e-^'(^^^-y^^t^'-l ^-^^ . (2.48) 

\ k — X k — ji J 

The analytic properties of this ■0 suggested the special choice for the Spectral Transform 
introduced in section 3. 



2.3 Wave soliton solutions of the Hamiltonian DSI 

To obtain soliton solutions of the Hamiltonian DSI equation ( |l.4| ) we use the elementary 
Backlund gauge in which we keep the general space-dependence of the parameters as 
indicated in ( 2.16 ) and ( ^.17 ). We start from the solution Q = as before but, now, we 
choose an auxiliary field A of the form diag(AgQ'' (u, i), Aqq"* (u, t)) where Aqq'' and A^q^ are 
two arbitrary functions. 

The construction of the soliton solution proceeds in the same way as previously, i.e. 
by imposing the commutation relation (2.22). By means of the elementary Backlund 
gauge of the first kind Bj{Qj, 0; A) defined in (2.16) we obtain the following solution of 
(1.5) (which does not satisfy the reduction) 

_fA^Piv,t) 



Qi = 





ri 



Ar = 









Af\u,t) 



where 



riiu,v,t) = Piu,t)/D{v,t), 



1(1), 
1(2), 



AV'{v,t) 



AV'{u,t) 



Ail\v,t) + 2dtlogD{v,t), 



(2.49) 

(2.50) 
(2.51) 
(2.52) 



The functions D(v,t) and P{u,t) are solutions of the time dependent Schrodinger equa- 
tions 



4(2) p _ 



0, 




(2.53) 
(2.54) 



and 



X{v,t) ^ d^\ogD{v,t). (2.55) 

By applying now to the above solution the elementary Backlund gauge of the second 
kind Bjj(Q, Qj; /i) we obtain the solution 

'0 



where 



Q 

q = 
A(2) 



A 



AW 








ED + PH/A ' 



= A 
= A, 



(1) 
00 

(2) 
00 



2dl\og{ED 
2dl \og{ED 



PEu - PuE 
ED + PH/A' 

PH/A), 

PH/A). 



(2.56) 



(2.57) 

(2.58) 
(2.59) 
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The functions H{v,t) and E{u,t) satisfy the time dependent Schrodinger equations 

ii7t+i7™ + 4o'^ = 0, (2.60) 
iEt-Euu + A'-;^E = {). (2.61) 

Note that the function ^ and the operator I = \ diag(9^ ^, d^^) are not uniquely deter- 
mined but are related by the equation 

^x{u,t) = -\d-\rq) - du\og{ED + PH/A). (2.62) 

The boundary values Ag^^ and Ag^^ defined in (1.4) are computed explicitly from the 
solutions of the four Schrodinger equations and from their potentials (the inputs A^^^ 
and 4'o)) 

= 4o + dl \og{ED + PH/A)\u=-oo + dl \og{ED + PH/A)\^=+oo, (2.63) 
4" - 4o - € \og{ED + PH/4)|,=_oo - dl \og{ED + PH/A)\,=+^. (2.64) 

One could get the same solution Q, A by applying the Backlund gauge of first and 
second kind in the reversed order Bj{Q, Qu; X)Bjj{Qii ,0; /x). One gets 

where 



and 



qii{u, V, t) = H{v, t)/E{u, t), (2.66) 

A^lj\v,t) ^ A^^„\v,t), (2.67) 

4^; («, t) = 4;;' («, - 29^ log e{u, t) (2.68) 

/i(u,t) = -a„logi;, (2.69) 



while A is not uniquely determined. If one requires that the commutation condition (2.22) 
is satisfied one gets the same solution Q^A, both A and are determined according to 
the equations ( 2.69 ), ( 2.55 ) and the operator X satisfies the equations 

du^qr = qri, d:^^ qr = -qur. (2.70) 

Let us consider the special case of arbitrary boundary values Agg real and moving with 
constant speed 

4;)(«,t)=A«(t; + 2^i) 

A'i^{u,t)^A'i^{u + 2et). (2.71) 

If the phases of the complex functions D, H and E, P are chosen to be linear functions 
it results that 

D = I'exp[-z(5], = ?7exp[-i(5], (2 72) 

E = £eiqy[ie\, P = pcxp[ze] 
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where 



and 



V = V{v + 24>t), r] = Tj{v + 2,^t), 

£ = £{u + 26t), p = p{u + 2et) ^ ' 



S = (j)v + {(j)^ - (t)o)t + 5o, 

e = eu + {e^- 9o)t + CO (2.74) 

with (1)0, $0, 6o, eo real constants. The real functions V, r] and £, p satisfy the Schrodinger 
equations 



V,, + {A'^^J -ct>o)V = 0, 
Vvv + (^00 - Mv = 

and 

£uu — (^00 + ^o)£ = 0, 
pun - (4o^ + Oo)p = 0. 

The solution Q is given by the following formulae 



(2.75) 

(2.76) 



where W is the wronskian operator. The field A is given by the formulae 

^(1) = A« + 292 iog{£V + pV4), 

A(^)=A^^^-2dl\og{£V + p7i/4) ^ • ' 



and its boundary values by 

4'^ = 4o + logiSV + pv/^)\u=-oo + dl \og{£V + pr;/4)|„=+oo, 
4'' = 4^0^ - dl \og[eV + p^|^%=-oo - dl \og{£V + prilA%=+^. 

The reduced case r = aoq is simply obtained by requiring that, for a G M, 



(2.79) 



W{r),V)^2a (2.80) 
W{p,£) = -2aoa. (2.81) 

In conclusion we obtain, in the reduced case, a solution of the DSI equation depending 
on two arbitrary real functions A'^J and A'^q . 

The one soliton solution exponentially decaying in the plane can be obtained from 
this general solution by choosing 

^ttd = ^Im =0' '^o = <^o, P-a = ^0, 

2? = exp[Ao(w + 20t)], = /iq cosh[Ao(w + 2(/)i)], (2.82) 

£ = exp[—pQ{u + 29t)], P = Po cosh[/xo(w + 20t)] 

where Aq, po, ho, po are real constants satisfying the constraints aoPoPo = Aq/iq and 
hoPo > 0. However, we are here interested in getting solutions with identically zero 
boundary values 

= A'i'> = (2.83) 
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in the reduced case r = dog. Because this case has been shown to be integrable in the 
Hamiltonian sense we call it, for brevity, the Hamiltonian case. 

It is convenient to introduce two functions 770 (w + 2(j)t) and po{u + 29t) defined as 
follows 



id) 



292logr,o, A, 



(2) 
00 



-25„logpo. 



(2.84) 



Then, to get a solution of the Hamiltonian case we have to solve the complicated non 
linear system of coupled equations for I?, ry, 770 and p, pq 



dl log r]oV 



AEV 
PV 



92 log 



u—-\-oo 



v—~{-oc 



= 0, (2.85) 
0. (2.86) 



Therefore we have to solve the non linear system of coupled equations ( 2.75| ), ( 2.76| ), 
(2.85) and (2.86) with the constraints (2.80) and ( |2.8l| ). We add to these equations the 
additional constraints 



lim ^ — 

u — »±oo c 

lim ^ = 



-2(pi±p2), 

-2(771 ±772) 



(2.87) 
(2.88) 



where pi and rji are real constants to be determined. These requirements allow us to 
decouple equations (HH), ( |2.75| ) and ( |2.80D from equations ( |2.86D , (|2.76D and ( |2.8l| ). 
Once having found a special solution P, 7/, 779 of the first group of equations and a 
solution 5, p, pq of the second group one has to verify that they satisfy, respectively, the 

requirements ( 2.87 ) and ( 2.8§| ). 

Let us first consider the equations ( |2.85D , ( p75| ), ( |2.80| ) and ( |2.87]). F rom ( ^.85] ) and 
(2.87) we get, by using the indeterminacy in the definition of 770 in ( ^.84 ), 



(2.89) 



If we express T) in terms of a new function a{v + (pt) as follows 



1 ,1 

2Pi'n + ^ 2P2T1 coth a. 



/2 , 



(2.90) 



we obtain 



2 sinh a 



V = 



P2770 

pi sinh a + tT'y02 cosh a 

P2??0 



where a and 770 are to be determined by requiring that (2.75) and (2.80) are satisfied. It 
results that 770 decouples from a 



(2.91) 
(2.92) 



dtm + a P2% - '/'oTyo = 



(2.93) 
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and a can be determined in terms of rjo 

dyU = a'ap2riQ. 



(2.94) 



For solving the equations (|2.86| ), ( 2.76 ) and ( 2.81 ) wc introduce, in an analogous way, a 
new function /3{u + 29t) as follows 



a"' = 1. 



We get 



P - 
£ 



2sinli/3 

rji sinh (3 + (j"rj2 cosh j3 
V2P0 



where po and /3 are determined by the equations 

dlpo + a^ripl - 6I0P0 = 0, 
duP = -aoa"ari2pl. 

Finally one has to verify that the following consistency conditions are satisfied 

lim (pi +(T'p2C0tha) = -(r^i ±7/2)""^, 

V — >±oo 

lim (771 + o-"?72 coth/3) = -(pi ± p2)"^ 

u — *±oo 

The solution q can be written as 

2ap2i]2PoVo exp[-i(e + 5)] 



q = 



sinh a sinh /3 + {pi sinh a + a' p2 cosh a) (771 sinh /3 + cr"?72 cosh 0) 



(2.95) 

(2.96) 
(2.97) 



(2.98) 
(2.99) 



(2.100) 
(2.101) 

(2.102) 



The ordinary differential equations ( ^.93 ) and ( ^.98 ) for rjo and po can be explicitly 
integrated in terms of elementary or classical tr ansce ndenta l func tions and, consequently, 
it is easy to verify the consistency conditions ( ^.100 ) and ( 2.101 ). 
For the sake of definiteness we consider two cases 



(i) dyr]Q = 0, 

(ii) d^rio ^ 0, 



In the case (i) 



duPo = 0, 
duPo ^ 0, 



2 '^o 
ap2 



(j)Q^\l> 0, 
4)0 < 0, 



= Mo > 0; 
00 <o. 



pI 



am 



and the consistency conditions are satisfied for 

a'a" XqHq > 0, pi± a' sgn(Ao)p2 = - 
If we choose for instance cr"Ao > we get the infinite wave 

q{u,v,t) = 2|AoPo| 



1 



m ± V2 



i/2 exp[-7(e + ^)] 



cosh^ 



(2.103) 



(2.104) 



(2.105) 



(2.106) 
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where 



^ = fio{u + 29t — uo) — Xo{v + 2(f)t — vo), 



In case (ii) equation ( 2.93 ) can be integrated once to the equation 



(2.107) 



(2.108) 



with 000 an arbitrary constant that we choose less than zero. Its general solution rio{v + 
2(f)t) can be expressed in terms of the Weierstrass elliptic function p{v + 20t; g2, 53) with 
invariants 

92 = ^1 53 = ^a'pl^lo " ^00 (2.109) 



and negative discriminant 

A(<72, 53) =51-2752 
according to the formula 



ia P20OO [a P2<l>t 



p{v + 20t;g2,53) - 0o/3' 



(2.110) 



(2.111) 



Note that 'i]o{v) for real v is always regular and that there exists a pure imaginary vq 
such that 

00 - . . 2 

3 ' 



p(wo;52,53) 



dvpi^o; 92,93) = i—^ap2(f>oo- 



(2.112) 



The function a{v + 2cj}t) obtained by integrating (2.94) (the invariants 52 and 53 are 
omitted and ao is a constant; C. and a are the (- and cr -Weierstrass functions) 



y{v + 2(j)t) = 



o-'ap20i 



'00 



dvpivo) 



2(x; + 0t)C(«o) + lof 



(^(w + 20t — Vq) 

a{v + 2(j}t + vo) 



+ ao 



results to be real and to behave for large u as follows 

a{v + 2(j)t) > 2a-'ap20ooTr^T^^^- 

Analogously we get for p^{u + 29t) {Oq < 0) 
Pq(u + 20t) = 



p{u + 2et;h2,h3) ^9o/3 
where the Weierstrass function p{u; h2, h^) has invariants 



(2.113) 



(2.114) 



(2.115) 



h3 = ^a^viOlo - 



(2.116) 
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negative discriminant 



A(/l2,/i3 



27hi 



and 



9 2 
P(mo; h2,h3) = —, dup{uo; h2, h^) = i--j=ari2dQo 



(2.117) 
(2.118) 



with uq pure imaginary. 

The function /3(u + 29t) is given by 



/3(u + 26li) 



a"aoar]20o 



2(u + 6li)C(Mo) + lof 



(t(u + 20t - uo) 



a{u + 2et + Uo) 



+ /3o (2.119) 



(invariants and /13 are omitted), which is real and has the following behaviour at large 
u 



/3(u + 2et) — > -2ao<7" ar]29, 



00 



u. 



dup{uo) 

The consistency conditions ( 2.10C| ) and ( |2.101 ) are satisfied for 



er' ct" p2doo4'oo 



and 



pi ± sgn ap20( 



dyp{vo)dupiuo) 
^ Civo) 



> 



1 



P2 



(2.120) 



(2.121) 



(2.122) 



By inserting these values and functions in the equation (2.102) for q we get an infinite 
wave with a periodically modulated amplitude. 



3 Solitons via the Spectral Transform 

3.1 The Spectral Transform of the Kadomtsev— Petviashvih I equa- 
tion 

It turns out that the Spectral Transform for the Kadomtsev-Petviashvili I equation plays 
a relevant role in the study of the DSI equation. Therefore, this section is dedicated to the 
main properties of this Spectral Transform that are of interest in this respect. Specifically, 
we will derive the multi-wave-soliton solution of the Kadomtsev-Petviashvili I equation 
and the orthogonality relations for the eigenfunctions of the associated spectral problem. 
We consider the Kadomtsev-Petviashvili equation in its variant (called KPI) 

{Ut - 6uUa: + Uxx:x)x = ^Uyy (3.1) 

with u = u{x, y, t) real. Its Spectral Transform is defined via the associated "time" (the 
space variable y plays here the role of time) dependent Schrodinger equation 

- + -u<P = 0. (3.2) 
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The spectral parameter k is introduced by requiring that 

$(a;,t/,fc)e''=^-*'='^ = l + O , A; ^ oo. (3.3) 

We consider, first, the case in which u is going sufficiently fast to zero at large distances 
in the (x, y) plane. Then, the eigenfunction $ can be chosen to be bounded in the {x, y) 
plane and scctionally holoniorphic in the complex k plane. More precisely. $ that is 
called the Jost solution is analytic in the upper and in the lower half plane and its 
boundary values on the two sides ±Imfc > of the real fc-axis are given by the 
integral equation 

^^x,y,k) = ldpe-'P-+'p"yR''{y,k,p), a = ±, (3.4) 



where 



R^iy, k,p) = 5{k-p)- sgn(p - k) 
X jdT]9{a{y-Tj){p-k)) jd^e'P^ 



'>uii,vWii,V,k). (3.5) 



When it is not differently indicated the integration is performed all along the real axis 
from — oo to +oo. 

The Spectral Transform F{k, I) of the potential u{x, y) is defined as the measure of 
the departure from analyticity of the Jost solution $ 

-= = / dl Adim)F{k,l) (3.6) 
dk J J 

where 

d_^lf_d_ . d \ 



- o I «; +*"^ ' fcffi = ReA;, fccj = Imfc (3.7) 
dk 2 \dksit dkcsj 

is the so called 9-derivative. 

The main quantity to study in order to get complete information on the Spectral 
Transform of the KPI equation is the scalar product 

($'^(fc),$^'(p)) = ^ jdx^%x,y,k)^-'{x,y,p). (3.8) 

By inserting here the integral equation for (f>'^ and <I>'^ we get 

($'^(fc),$-'(p)) = (R-(y)R-'^(y))(fc,p), (3.9) 

where R^^y, k,p) is considered as the kernel of an integral operator R'^(y), R'^\y, k,p) = 
R'^iViPi k) is the kernel of the adjoint operator R'^^(j/) and 

(R^y)R-'\y)){k,p)= jdl R%y,k,l)R-' {y,p,l) (3.10) 



20 



is the kernel of the product R'^fylR'^ (y). By differentiating the scalar product (3^ 



with respect to y and by using (3^) for and one proves that it is y independent. 
For exploiting this information it is convenient to introduce 



i?J(fc,p)= lim R^[y,k,p). 

y—^±oo 



(3.11) 

They are kernels of the triangular integral operators R^. whose explicit expressions 



r"(fc,p) = ^ Jd7jJd^e'Pi''P"'^u{tvW{^,V,k), a = ±, 

are obtained by inserting into (^^) the identity 

asgn{p-k)d{cr{y-ri){p-k)) =^d{T{y-ri)) ±d{Ta{k-p)). 
Then, the y independence of the scalar product implies that 

($-(fc),$-'(p)) = (R;R;'^)(fc,p) = (R^R^'^)(fc,p). 



(3.12) 
(3.13) 

(3.14) 

(3.15) 



In the case a' = —a the two operators RJ^R!^'^^ and R^RI'^^ are, respectively, lower 
and upper triangular or upper and lower triangular according to the sign, positive or 
negative, of a. Consequently, the equality in (|3.15 ) implies that 



R(T Tj — CT T 



R^RI'"^ = I, 



(3.16) 



where the kernel of the unity operator is I{k,p) = 5(k — p), so that and $ are 
orthogonal 

{^^k),^-^p)) = S{k~p). (3.17) 

In view of the relevant role played by R^ it is convenient, by using again the identity 
( 3.14 ) and the definition ( 3.11 ), to recast the integral equation (3^) in the form 

^"{x^y^k) ^ Jdpe-'P''+'P'yRl{k,p) 

-77- Apr d?7fee'f(«-^)-*p'(''~J^)u(e,?7)$'"(C,?y,fc). (3.18) 
27ri J J±oo J 

Let us now turn to the eigenfunctions solution of the integral equation 

^±{x,y,k) ^ e-'''''+''''y - (dpT dryfee'P(«-^)-*P'("-^)M(C,?7)*±(^,r;,/c). (3.19) 

27ri J J±oo J 



From (3.18) we see that 

= R±*±. (3.20) 

As we know by ( 3.16| ) the operators Rl^"^^ are right inverse of the RJ. Let the potential 
u be such that these operators are both sides mutually inverse, i.e. let be besides (3.16|) 



RT'^^Rl = R^^RI'" 



(3.21) 
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Then the relation ( 3.20| ) can be rewritten as 



(3.22) 



Due to the definin g int egral equation ( 3.19 ) the solutions ^>± are a independent. Conse- 
quently equation (3.22) yields 



and, by using again {3. It), 
where 



^—cr r>cr-r>crT 

J- _ rv_j_rv_j_ 



(3.23) 

(3.24) 
(3.25) 



From (3.15) and (3.25) we have directly that 

(<I)-(fc),$-(p)) ^^-(fc,p). (3.26) 
It is convenient to separate the 6 distribution contained in the kernel of IF'^ by writing 

T^{k,p)^6{k-p)-or{k,p). (3.27) 
Then, equations ( 3.24 ) and ( 3.25| ) read 

$+(x,2/,fc)-$-(a;,y,fc)= Jdl T {kj^ ix,y,l), (3.28) 



f-'^{k,p) = r''{k,p)d{k -p) + r'^ip, k)d{p -~k)+a dM{p- l^ik, l^^ip, I). (3.29) 

J —OG 

By recalling that on the real A:-axis 

^(fc)^^($+(fc)-$-(fc))<5(fca) (3.30) 



we get for the Spectral Transform of u defined in (3.6) 



F{k,i)^-r{kj)sik.^). 



(3.31) 



Formula (3.29) solves the direct spectral problem furnishing an explicit expression of 
the Spectral Transform in terms of u and of the Jost solution via r°'{k,p). 

We note that the operators JF'^((t = ±), which yield the Spectral Transform of u, are 
selfadjoint 

= (3.32) 



and one is the inverse of the other 



(T -T- — cr 



(3.33) 



These two equations can be considered as the characterization equations for the Spectral 
Transform !F" (in the sense that they determine the class of admissible spectral data) 
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if T'^ is the product, as in (3.16), of two operators which have the triangular form 
indicated in ( |3.12 ). Because the characterization equations are more simply expressed in 
terms of we prefer, in the following, to define as the Spectral Transform of u. 

The reconstruction of the potential u can be performed starting from the Spectral 
Data J-" . First, one solves the singular Fredholm equation 



^"{x^y.k) = e 



. —ikx-\-ik y 



dq ■ 



-i(k — q)x-\-i{k^ —q'^)y 



2Tii J ^ q — k — iOcr 

(3.34) 

The obtained by the Cauchy-Green theorem, solves the non loc al Ri emann-Hilbert 
problem expressed by the 9-equation (3.(;) (or its equivalent form (3.24)) and satisfies 
the asymptotic requirement ( [3.3[ ). Secondly, one expands in powers of 1/fc the fraction 
l/{q — k — iO(j) in ( 3.34 ) and inserts the obtained asymptotic expansion in (3^). We have 



u{x,y) = -a'^-^ jdqjdplJ'^' -l]{q,pW' {x,y,p)e^'i^-^^'y. (3.35) 

In the following section we are interested also in potentials u describing N interacting 
wave solitons. Then, u goes to a constant along N directions in the (x, y) plane and 
the Green function in the integral equation (3.4) has to be corrected in order to avoid 



divergences. The theory of the Spectral Transform does not result to be substantially 
different from the theory in the case of u vanishing at large distances. However, there 
are many subtle technical difficulties to handle which are out of the scope of the present 
paper. Therefore, we restrict the discussion to the main points and we refer for details 
to the published papers by two of the authors of this paper (M. B. and F. P.) and by 
A. Pogrebkov and M. Polivanov. 

The Jost solution v{x, y, k) — <E>(x, y, k) explikx — ik'^y] in the complex /c-plane can 
be defined as the solution of the following integral equation 



v{x,y,k) = 1 + 



1 



dpdq ■ 



1 



(27r)2 J J ^ q-{p- iO){p + 2k) 



(3.36) 



where the integrations must be done in the order indicated from the right to the left. For 
a potential u going sufficiently fast to zero for large x^ + y^ the order is unessential and 
one can explicitly perform, first, the integration over q recovering the integral equation 
0- 



The integral equation (3.36) contains as a special case the unidimcnsional case. In 
fact if u(x,y) = u{x) by introducing 



$(a;,fc) e-'''^y<S>{x,y,k) 
one recovers the stationary Schrodinger equation 

^xxix, k) + (fc^ - u{x))'^{x, k) = 0. 



(3.37) 



(3.38) 



Moreover, in (3.36), once integrated over 77, the integration over ^ and p, q can be inter- 
changed and one obtains by computing the limit fccj — s- on the two sides of the real 
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k-axis the familiar integral equations 



$ (x, k) ^ e 



— ikx 



0{i - xf^^^^^^-^u{0<^- k) (3.40) 



defining the sectionally meromorphic Jost solution $ of the stationary Schrodinger equa- 
tion. 

In the general bidimensional case, as in the previous case {u vanishing at infinity), the 
main quantities to study are the scalar pro ducts of the Jost solutions. They result to be 
still y independent and, therefore, formula ( 3.15 ) remains valid. The most difficult point 
is the co mput ation of the limits i?J and of the characterization equations. It results that 
formula ( 3.12 ) must be changed as follows 

RUk,p) = S{k-p){l + Zl{k)) T e{TcT{k-p))r%k,p) (3.41) 

by adding to the coefficient of the 6 distribution a function Z^{k) to be determined. 

The Jost solutions <I'°' and <I>~°' are still orthogonal and this is the main property we 
will use in the following section. 

Let us now proceed and consider the iV-soliton solution of the KPI, namely the 
potential u characterized by a Spectral Transform containing exclusively a number of 
discrete eigenvalues G C (n = 1, . . . , N). For real u this spectrum fd{k, I) must satisfy 
the characterization equation 

7,ik,l) = Ul,k) (3.42) 
and therefore its most general form is 

N N 

fd{k, l)=Y.Y. ^"™^(' - T^rn)5{k - fin) (3.43) 
n— 1 m— 1 

with r„m an arbitrary constant hermitian matrix. It is convenient to parameterize this 
matrix as follows 

^nm — 27rZ exp[i^^^(Xo77i ^'myoni) ^^n{,'^on f^nl/on^^Cnnip'rn'm (3.44) 

where the real parameters Xon and yon fix the initial position of the n*'* soliton in the 
plane, 

'Pnm = Mn ^ l^m, (3.45) 

and the complex matrix C satisfies 

^nn ~ ^nrnPmm ~ ^mn'Pnn- (3.46) 

For definiteness we choose C„„ = 1 and Im /i„ > 0. 
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Because the Jost solution $(a:, y, k) in the pure discrete case has only simple poles 
at k — fin [n — 1,...,N) and satisfies the asymptotic property (3.3) it admits the 
representation 

*(x, y, k) = e-^'=-+^'='« I 1 + E '^"^'"'^^ e"" ) (3.47) 



( ^ 
ikx+ik^ y I ]^ _|_ ^ ^ 



where, for convenience, it has been introduced an exponential factor e"" with 

Ctn = ifJ-nix - Xon) - ifJ-liy - Von) (3.48) 

The insertion of this representation into the 9— equation ( |3.6| ) yields an algebraic equation 
for the functions (pn 

N N 

E ^nmifrn = - ^ C'«™^mme"" (3.49) 
m— 1 m— 1 

where 

= 1 + Ca, a„™ = ft™ 5""+"™ , (3.50) 



while its insertion into (3.35) yields the potential 



u 



N 



= -2zc',^(p„e"". (3.51) 



To see that (3.49) can be solved for (p„ it is sufficient to show that A is positive definite. 
In the following section in the more general framework of the Spectral Transform for the 
DSI equation we prove that this is true if the hermitian matrix iCnrnptmrn positive 
eigenvalues. 

By the rule for differentiating a determinant the N soliton solution u can be written 
in the closed form 

~2dl\iidctA (3.52) 
and the orthogonality relations for the Jost solutions 



i dx^ (a;, y,/i„) Res (<l)+(a;, y, /c), = 5,„„ (3.53) 



can be derived. 

The potential u{x, y) is going to a constant at large distance along the directions 
X — 2Re /i„y = const and describes N intersecting waves of infinite length. 



3.2 The Spectral Transform of the Davey-Stewartson I equation 

We consider the DSI equation (|l.5| ) with boundary conditions defined as in ( |l.8| ). The 
real boundaries a^^^v^t) and a^\u,t) are assumed to go to zero at large distances in 
the {v,t) and (ii,t) plane, respectively, with the possible exception of a finite number of 
directions along which they are going to some constants. 

According to the usual scheme, in order to linearize the DSI equation, we have to 
define the Spectral Transform for the Zakharov-Shabat spectral problem in the plane 
(hyperbolic case) 

Ti^ = {d, + a:idy + Q)^/' = (3.54) 
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The complex spectral parameter k is introduced by requiring that the 2x2 matrix Jost 
solution tp satisfies the asymptotic property 

Mx,y,k)e-''''^'''''-y^ = l + 0(h, k ^ oo (3.55) 

k 

The Green function of the Zakharov-Shabat spectral problem can be chosen to be sec- 
tionally holomorphic and the values V'* of ip on the two sides ±Im fc > of the real axis 
in the fc-plane are given by the integral equations 



where 



= + G±V^ (3.56) 
G±V= (Gft/.i,G±t/>2) (3.57) 



4' = {tPl,tp2) 



^11 V'12 
V'21 '022 



with 



1 J^°°du'q{u',v) 

2 i;^°°dv'r{u,v') 

l( J^'^du'qiu'^v) 



2 V 1 dv'r{u,v') J 

and ip^ — diag(V'^]^, "002) arbitrary solution of the homogeneous part of the Zakharov- 
Shabat spectral equation 

(a, + asdy)^'^ - (3.59) 
It is worth noting for future use that the Green operators G^ have the symmetric property 

G+ = Gf (3.60) 

and are k independent. 

The 2x2 matrix Spectral Transform R{k, I)) of Q{x, y) is defined as the measure of 
the departure from analyticity of ip 

^tP{x, y, k) = jjdl A d7 V'(x, V, l)R{k, I). (3.61) 

In contrast with the one dimensional case, once chosen the Green operator, the asymp- 
totic requirement ( 3.55|) d oes not fix 0oi and which are arbitrary functions of v and 
It, respectively (see ( 3.59| )). The refor e, for different choices of 0o wc get different ip and 
consequently via the definition ( ^.61 ) different Spectral Transforms R{k,l). 

We search for a Spectral Transform that satisfies the following two fundamental prop- 
erties: 

i) its time evolution can be explicitly integrated; 

ii) the discrete part of the spectrum corresponds to solitons and the continuous part 
to radiation. 
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In order to satisfy the requirement i), in analogy with the one dimensional case, we fix 
the time evolution of i/' by requiring that 



(3.62) 



where T2 is the second Lax operator in (1.13) and fl{k) ~ —a^k^ is the dispersive function 
of the DSI equation. In fact, by applying to both sides of (3.62) the operator ^ and by 
using the definition of the Spectral Transform in (3.61) we get the linear time evolution 
equation for the Spectral Transform 



iRt{k, I, t) = R{k, I, t)n{k) - n{l)R{k, I, t) 
which is easily explicitly solved getting 

R{k, I, t) = e^"(')*i?(fc, /, 0)e-*"W* 



(3.63) 



(3.64) 



By taking the limit (u ~oo, v fixed) for the ( )ii matrix element of (3.62) and the 
limit (w —^ —00, u fixed) for the ( )22 element we derive two time dependent Schrodinger 
equations for ^/jqi and ipQ2 



[dl + fc2 + a\^^ (v, t)] Vol {v, t, k) = -ift Voi(v, t, k) 
[dl + fc2 - a^;^\u,t)]^oi{u,t:k) = idti>02{u,t,k) 

These equations together with the asymptotic property 



]i + o(-), 



00, 



(3.65) 



(3.66) 



which can be derived from (3.56) by using (3.55), uniquely determine f/'o in terms of the 

Precisely, we have 



boundary values Oq^^* and Cq 



V'oi (i^,i, A;) = $(!)(«, t,fc)e^'='* 



(3.67) 
(3.68) 



where (w, t, k) and <I>^^^(u, t, k) are the Jost solutions of the time dependent Schrodinger 
equations with space varia bles v and u and potentials a^^^ {v, t) and a*^^-* {u, t), respectively. 

The integral equations ( 3.56| ) for ip are of Volterra type and therefore the singularities 
of V' in the complex /c-plane are those of V'o and of the sectionally holomorphic Green 
function. Therefore we need to revisit the Spectral Transform for the time dependent 
Schrodinger equation. In particular if ag^'' and Oq^^ are wave solitons in the plane (v,t) 
and (m, t) the eigenfunction tjjQ and consequently ip have simple poles in k. The existence 
and the location of poles are uniquely determined by the boundary values Oq^^ and a^'^K 
However, we will see that the boundaries do not characterize completely the discrete 
spectrum and one is left with an unexpected freedom in choosing other independent 
parameters. 

The continuous part Rc{k,l) of the Spectral Transform R{k,l) measures the discon- 
tinuity -0+ — ?/;_ of ^ along the real k axis and it has the form 



Rc{kJ) = — 



1 f S{h + 0) 






6ih - 0) 



FiikJ) 
Si{k,l) 



-S2{k,l) 

F2{k,l) 



S{k. 



(3.69) 
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where we make the choice ( "^^^^ + ° ^ instead of the usual one S(L> + 0)1 

\^ o[Iq-0) J 

or S{Iq — 0)1 in order to exploit in the following the symmetry ~ of the Green 

operator. 

Let us introduce the integral operator 

G= (G^,G+). (3.70) 

From the integral equation ( 3.56| ) we derive directly (fc is real, space time variables are 
understood) 

(^+ _ = j^^ik) + G(^+ - V~)(fc) (3.71) 

where 

^c(fc) = . (3.72) 

\ ^Jdvri;+,ik) (^o+2-V^2)(fc) / 

Here and in the following when it is not differently indicated the integration is performed 
all along the real axis from — oo to +oo. 

By inserting the same integral equation (3.56) into the r.h.s. of (3.61), by using the 
symmetry G2 = and the k independence of G* and by recalling that on the real k 
axis 

^(fc) = l(^+(fc)-^-(fc))5(fccj) (3.73) 
ok i 

we get ^ 

(^+ _ ^ + G'(^+ - i,~){k) (3.74) 

where 

J',{k)= dl\ (3.75) 

J \ ^o+(05i(fc,0 ^0^2(0^2(^,0 J 

Since the operator G is of Volterra type the related homogeneous integral equation has 
only the vanishing solution and, consequently, by comparing the two integral equations 
(|7l]) and (|j|) we get 

= (3.76) 
From this equation we deduce that Fi and F2 are the continuous component of the Spec- 
tral Transform of the boundaries a^Q^^ and Oq^'' and, by using the orthogonality relations 
(3.17) rewritten for the eigenfunctions ipoi and ipQ2j we express explicitly 5*1 and 52 in 
terms oi Q, tp and ipQ as follows 

Sl{k, — ~ J J dv r{u, u)-0j*'i("j ^: ^)V''02('"' 

■^2(^,0 -^j jdudvq{u,v)il;22{u,v,k)^^^{v,l) (3.77) 

In order to get the characterization equation satisfied by the spectral data Si{k, I) in the 
reduced case r = a^q we consider the two Jost solutions 



^2ik)=( ^-^(fc) ]+G2rP2{k) (3.78) 
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and we note that from (3.54), for r = aoq, it easily follows that 

(V'+i (O^(fc))^ = CIO (^2+ (OV^(fc))^ . 



By integrating it in the (u,w) plane and by using ( 3.78 ) we get 



S2ik,l) crQSi{l,k). 



From (3.64) we derive that the spectral data evolve in time as follows 

Siik,l,t) = e*('='+'')*S'i(fc,/,0) 



(3.79) 



(3.80) 



(3.81) 
(3.82) 



coherently with the characterization equation ( 3.80| ). 

The discrete component Rd{k, I) of the Spectral Transform has different possible 
characterization. We consider two of them whose matrix elements are linear combinations 
of S distributions in the complex fc-plane and /-plane. 

The first one is given by the formula 



Rd{k,l) = -27ri 

N N 



Rii{k,l) i?i2(fc,0 

R2l{k,l) i?22(fc,0 



(3.83) 



Ruik,l) = ^ ^ T^r'Sil - ^n')5{k ~ A„) 

n— 1 n' — 1 

Rl2{k,l) = Y.T. - ^n)S{k - /i„,) 

n—1 m—1 
N M 

R2i{k, = E E ^2T^a - J^rn)S{k - A„) 



n—1 m—1 
M M 



R22{k,l)^J2 E 



'22 



^{1 - tJ-m')Hk ~ Aim) 



m — 1 m' — l 



where A„ and /i^, the so called discrete values of the spectrum, are, respectively, the 
locations of poles in the complex fc-plane of the first and second column of the Jost 
matrix solution and the t^™" are some complex constants to be related to the initial 
value of Q. 

The S distribution in the complex plane is defined as 



dlAdlS{l-lo)f{l) = f{lo) 



(3.84) 



If we introduce the distribution (I — lo)S{l — lo) operating on singular functions in the 
complex plane as follows 



dl Adl{l~ Iq)S{1 - lo).f{l) = Res(/, ^o) 



(3.85) 
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we can deal with a simpler form of the Spectral Transform characterized by an off diagonal 
matrix 

N M 

n—l m— 1 
Af M 

R2{k,l) = Y.Y1 [rni - >^n)S{l - Xn)S{k - fi^) + ?rd{l -K)6{k - ^i^)] . 



By inserting the two different characterization of Rd in the i9-equation ( |3.61| ) one derives 
easily that they are equivalent and, more precisely, that the two set of spectral data are 
related by the following equations 

Til = -(4 + rir2)"Vif2 
Ti2 = —(4 + r2ri)^"'2f2 
T21 = (4 + riT2)"^2n 

T22 = -(4 + r2ri)-V2n (3.87) 

and 

Tl = 2tiiT^2 

T2 = -'2.T22T21 

Tl = 2t21 — 2tiiTi2^T22 

T2 = -2ti2 + 2X22X21 Vll. (3.88) 

Without any loss of generality we can consider M ^ N . The cases M < N and M > N 
are recovered from the special case M = A'^ by choosing some /i or A equal. 

In order to distinguish between parameters that are fixed by the choice of the bound- 
aries, i.e. the external data given at all times, and parameters that are solely connected 
t o th e initial data at time < = it is convenient to parameterize the spectral data in 
( 3.86 ) as follows 

t"" = Pnm exp[i^„iUo„i + iXnVon + i^J-l^t + iX^t] 
N 

~nm ^ ^ PnpCpm-Pmm exp\ijL^Uom + iXnVon + ^Mm^ + 
p=l 

T^" = Vrnn exp[-«^„iU<„„ - iX„Von " ^Mm* " «^«*] 
N 

~mn ^ ^ rjjnpDpn-Xnn exp[-l^„Uo„ - iXnVon ~ ip-^t - iX^t] (3.89) 
p=l 

where the time dependence is explicitly given, 

'Pmn ~ f^n 

-Xmn = Xm ~ A„, (3.90) 
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p and 77 are arbitrary complex constant matrices, Von and Uon are arbitrary real constants. 
We shall show that the choice of the boundaries Og^^ and Oq^^ determines uniquely the 
complex matrices C and D, together with the A and the /i, and that the other parameters 
are left free. For definiteness we choose 

ImA„<0, lm^„>0. (3.91) 

By computing the residua at the poles fc = A„ and k — /i„ of both sides of the integral 
equations ( 3.56| ) we get 

Res(V^i,A„) = I^es(V'oi,A„) j +GrRes(V^i,A„) 

Res(^2,Mn) = ( Res(^lA.„) ) + ^2+ Rcs(^2, (3.92) 



By inserting these integral equations into the r.h.s. of equation ( 3.61 ) considered at the 
special values k — Xn and k — /i„ and by recalling that the 9-derivative of a pole at 
A; = fco is given by 

^I-^ ^-27rz,5(fc-fco) (3.93) 
dkk-kQ 

and the symmetry property of the Green operator G 

G2 = 

_ 1/ Jduq{u,v) 

^2 +^[0 



we get, for the special choice of Rd{k, I) in ( |3.86 ), 



Res('0i, A„) = J^i+G^ Res{ipi,Xn) 

Res('i/'2:Mn) = ^2 + GjRes(i/'2>Mn) (3.95) 



where 
with 

AT 



^1 = ( ?M (3.96) 



•^11 = ^ X! A^'?V'22(Mm)(pC')n™t*™mexp[i/x„Wo^ + iA„i;o„ + ZM^i + iA^i] 

m=r 

1 ^ 

-^21 2 X! {F'''5s(V'02; Mm) exp[iyUmMom + iXnVon + il4nt + + 

m— 1 

CmpPpp1p02 i'Pp) exp[iJIpUop + iXnVon + ip^pt + iX^t] } (3.97) 



0=1 
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and 



Ti = 



12 



22 



(3.98) 



with 



1 ^ 

^12 = - X ^ ?7nm {Res(V'oi, Am) exp[-i/i„Mo„ - lAmUom - ~ «^m^] + 



TV 



p=l 



1} 



(3.99) 



By comparing the integral equations (3.92) and (3.95) one gets 
Res(V',7i, A„) = 

- ^ /(iugV22(Aim)(pC')nmt*mmexp[i/X 



(3.100) 



m— r 

',+ 



Res(V'o2,Ai«) 



- ^ I dv r^^^{\m){'nD)nnAmmey^v[-i 



^^JinUon - iXmVom - ^^J.nt " ^KA (3.101) 



and 



N 



Res(V'oi,A„) = - ^ Dnrn^rnmlpoiiXm) exp[iX iX rn'^om 

+ iXlt^iX„A (3.102) 



N 



Res(V'027 Mn) = X! <^™"^mm^02(Aim) exp[-i^„Uon + ifJ-m'^o^n " iMn* + ^A^m^] (3.103) 
m— 1 

Equations (3.102) and ( 3.103| ) determine uniquely the boundary values and a^^^ of 
the auxiliary field A = dia,g{A'^^\ A^^''). In fact, they determine the Spectral Transform 
of a^^ ^ and a'i^^ considered as potentials in the time dependent Schrodinger equations of 
(3.65). If we use the Spectral Transform defined in the previous section 



dtpo 



dk 
dk 



dl /\dli:oi{l)ri{k,l)e''-'' 
Ad7Vo2(0^2(fc,0e~*''''"''^ 



(3.104) 
(3.105) 



we get that 

ri(fc, I) = 27ri ^ exp[iA„Uo„ - iX„iVom]DnrrAmmS{l - X„i)S{k - A„) (3.106) 



N 



n,m— 1 
N 



r2 



(fc, I) = 27ri ^ exp[-i^„Uo„ + l^„Uo„]C„„t«mm'5(^ " f^m)^(k - fJ-n) (3.107) 
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It results, in particular, that i/'oi smd 7/^02 have simple poles at fc = A„ and at k = fXn as 
expected. 

If we consider the KPI equation associated, for instance, to the time dependent 
Schrodinger equation for a'l^\v,t) by taking into account that {v,t) are to be considered 
as "space" variables of the KPI equation while its "time" variable has to be considered as 
an additional parameter of the potential al^\v,t), it results that the obtained boundary 

value a'j^\v,t) coincides with an N wave soliton solution of the KPI equation in the 

(2) 

{v,t) plane at some fixed "time". Analogously, Oq {u,t) can be considered as an N wave 
soliton solution of a KPI equation in the (u, t) "plane" at some fixed "time" . The real 

th 

parameters Von and Uon fix the position of the n wave soliton in the corresponding 
"plane". 

If we require Oq^'' and ag^^ to be real the complex matrices C and D satisfy 

Dnrn^mra — Drart^nn (3.108) 
Cnmpmm ~ ^mn-pnu- (3.109) 

Moreover, without any loss of generality, by eventually shifting Uon and Von we can choose 

dL = CL - 1 (3.110) 



Equations ( 3. IOC ) and ( ^3.101 ) solve the direct proble m fur nishing the matrices p and 77 
in terms of the eigenfunction of the spectral problem ( ^.54| ), of Q and of the parameters 
and eigenfunctions defining the boundary. By using the orthogonality relations derived 
in the previous section 



j/dt; V'oi(Am)Res(-0oi' = ^r, 



«/t^U^02(/"m)R-es(7/;(j2,M«) =Sr, 



(3.111) 



we obtain the formulae 



Pnm = -iexp[-iXnVon " «A^t] 

Vnm = -?exp[i^„Uo„ + ifj.lt](Y^'^) 



(3.112) 



where 



N 



nrriPmrn^ 



ni—1 
N 



ciudi;rV'n(A„)V'o2(Ms 



(3.113) 



It is worth noting that the chosen parameterization of the spectral data in (3.89) allows 



us to discriminate between the parameters A„, /x„, Wo„, Vom Cmm Dmn which are fixed 
by the choice of the boundary values o'q^ and Oq^"* and parameters pmn, Vmn which are 
left free and are expected to govern the specific nonlinear dynamic of the DSI equation. 
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In analogy with the one dimensional case we call pmn and rjmn the normalization matrix 
coefficients. 

In the reduced case r — aoq from equation (3.8C) computed &t k = fin and at / = Am, 
by using the orthogonality relations (3.111), we get the necessary condition 

N N 

^ ] PnmCmsPss ~ '^0 ^ ^ D nrri^vamfl sm (3.114) 
m—1 m—1 

In solving the inverse problem we shall prove that this condition is also sufficient for 
having r = aoq. 

We are left with the solution of the inverse problem, i.e. we have to reduce the 
reconstruction of the matrix Q and the auxiliary field A from given spectral data i?(fc, I) 
to the solution of a linear problem. 

The 9-equation (3.61) together with the asymptotic requirement in ( 3.55| ) defines a 
non-local Riemann-Hilbert problem for the matrix function 

<l>{x, y, t, k) - ^(x, y, t, k)e~'''^^'^-y'> (3.115) 

Its solution is obtained by solving the singular linear integral equation (the space-time 
variables of are understood) 

</.(fc) - 1 + // ^ ffdlA d7</)(0e''("^"-^)i?(/i, I, t)e-^''("^"-'') (3.116) 
27r«J h ^ k J 

This equation furnishes an asymptotic expansion in powers of -i of 



(3.117) 



dhAdh dlAdl 0(Oe*'('"^'^-^)i?(/i, I, i)e-'''(''^^-«). (3.118) 



where 

Then the solution of the inverse problem is achieved by inserting the expansion into 
the Zakharov-Shabat spectral problem Ti'0 = and into the auxiliary spectral problem 
T2'0 = and by identifying the coefficients of the powers of ^. One obtains Q, A and 
a useful expression for 

g = ^[(73,(/.(l)] (3.119) 
A = -i(9,-a35^)diag0W (3.120) 
g2 = 2ia3idx +cr3dy)diag(t)'^^\ (3.121) 

When only the discrete part of the spectrum is present one can derive explicit algebraic 
formulae for Q and A. The requirement that (j) has simple poles at fc = A„ and at fc = /i„ 
and that (j) goes to II in the large k limit fixes a k dependence of the form 
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where, for convenience, in the residua an exphcit exponential factor has been introduced 
with 



an = ifJ-niu - Uon) " IMn* 
(in = -i\n{v - Von) + i>^lt 



(3.123) 
(3.124) 



The vectors ipni a nd are com puted by inserting ( 3. 122 ) into the (9-equation ( 3.61 ) 
with R{k,l) from (3.86) and (3.89). One obtains 



l+-vBpA 
l+lpAr,B 



= -7;vs(l] - -,vBpi 



2' \0J 4 



1 /0\ 1 



2'^'Vl/ 4 



-pAr]S 



where 



AT 

m— 1 
N _ 

m— 1 

^ = 1 + Ca, Q;„m 



From equation ( 3.119 ) one gets 



N 



Q = 2ia3 J2 



(¥'1)26"" 
(¥'2)1 e'^" 



(3.125) 
(3.126) 

(3.127) 

(3.128) 
(3.129) 
(3.130) 

(3.131) 



and from ( |3.120 ) and ( 3.121 ) by the rule for differentiating a determinant one gets 







)2 = -49,9„lnA 



In A 



where 



A = dot 



l+-vBpA 



= dot 



t+-pAvB 



(3.132) 
(3.133) 

(3.134) 



By solving, in a similar way, the 9-equations ( |3.104D and ( ^.105] ) one can derive also the 
boundaries 



= 29^1ndetS 



,(2) 



-2a,lndet A 



(3.135) 
(3.136) 
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In order to study the reduction and the regularity properties of the found sohition it is 
convenient to introduce the hermitian matrices 



One can easily get that 



CA^ = AC 
DB^ = BD 



Moreover, if the constraint (3.114), which rewritten by using C and D reads 

Dri^ = -cTopC, 

is satisfied one can easily verify directly that the reduction condition 

r = aoq 



(3.137) 
(3.138) 

(3.139) 
(3.140) 

(3.141) 
(3.142) 



is fulfilled. We conclude that the constraint (3.114) is necessary and sufficient for having 
the reduction. 

To see that ( |3l25| ) and ( |3l2^ ) can be solved for ipi and and that, consequently, 
Q and A are regular, we first note that the two hermitian matrices 



an 



-i 



(3.143) 



(3.144) 



are positive definite. The corresponding quadratic forms in the dummy vector P can be 
written as 



J2 PnPmanm^ ^"^™/ ^w' e""+"" = / du' 

N N 



N 



n,m— 1 



n,m — l 



71=1 



(3.145) 



(3.146) 



which are positive unless all P„ are zero. Therefore if the matrices C and D are chosen 
to be positive definite, because the product of two positive hermitian matrices is positive 
A = 1 + Ca — t + Ca and B = l + Df] = 1 + D(3 are positive definite and the boundaries 
Cq^^ and g^^^ ar e regu lar. I f, moreover, the reduction condition ( 3.137 ) is satisfied the 



matrices ( |3.125 ) and (3.126) can be rewritten as 

1 



^riBpA = 1 - ^aQr^BDr]'^C-^A 



1 + ^pAr]B = 11 - ^aapACp'^D-^B 



(3.147) 
(3.148) 



which for = —\ are positive definite for regular 77 and p. For ctq = 1 they are positive 
definite for rj and p matrices with sufficiently small norms. 
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3.3 The dynamics of localized solitons in two dimensions 

The solution q, in the generic case, describes N'^ locahzed coherent structures interacting 
in a compUcated way at finite times, but moving in the far past and in the far future with 
constant velocities V„„i = (2Re A„, 2Re /i„i) and without changing form. It is therefore 
natural to call this solution the N'^ soliton solution. It is parameterized by a point in 
a space of 4:N{N + 1) real parameters. Of these parameters 2N{N + 2) and, precisely, 
A„, i^n, Von, C and D are determined by the choice of the boundaries and fix the 
velocity and the possible location of the solitons in the plane, while the remaining 27V^ 
and, precisely, and p govern the dynamics of the solitons during the interaction. 

A relevant information in the global dynamical behaviour of the solitons is furnished 
by the mass (energy, charge, or number of particles according to the physical context) of 
the solution q 

M=ff\q\^dudv (3.149) 



and by the masses of the solitons at t — ±oo 

M^r! =ll\ <lin2 I' dudv (m, n = 1, 2, . . . , iV) (3.150) 



where qmn is the asymptotic behaviour of q at t ~ ±oo computed in the rest reference 
frame of the (m, n) soliton. If dct(p77) ^ 0, dot a and det /3 ^ the mass AI is given 
by the simple formula 

det (1+ hip) 

M = -4ao\n (3.151) 
det [jr]p) 

In general, it results that M = J2mn =J2mn ^^mn- However, the mass of the single 

soliton is not conserved and in particular it can be zero at t = +oo or at i = — cx). For 
a special choice of the parameters the mass of the (to, n) soliton can be zero at t = ±oo 
also when the coefficients and rjnm are not both equal to zero. We call these solitons 
with zero mass virtual solitons and they generate peculiar effects as in figure 1 of |^ 
where a virtual soliton collides with a soliton forcing it to change velocity. 
On the other hand, the total momentum of q 

P = (Pi,P2), Pi = iJJ {q.^q -qqu)dudv, P2 = iJJ {q^q ~ qqv)du dv (3.152) 

is not conserved and, in fact, it results that 

^^(^[f \q\^ iA2iu,t)dudv,-2 [[ \q\^ -^Ai{v,t)dudv) (3.153) 

where 



dt \ J J du ' ' ././ dv 



Ai = ai'^ - i^o / du{\ q p)^, A2 = + / dv{\ q p)^. (3.154) 

Because the four-soliton solution displays all the richness of the general case we 
examine it in detail. We choose for definiteness A„3, < 0, /i„Q. > (n ~ 1,2), (A2SR — 
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Aisr) > and (/i2K — Misr) < 0- It is convenient to parameterize the matrices p and 77 in 
such a way that the reduction condition (3.114) is automatically satisfied, i.e. 



V ^11 V ^22 



c 



-^22 I 



V 1^11 V ^22 



V 



V ^22 / 



(3.155) 



(3.156) 



with In and to„ (n — 1,2) arbitrary complex parameters and 



C = C 



12 




V = D 



1- 



f'li 



12 




1- 



^11 



(3.157) 
(3.158) 

The parameters A„ , /i„ , C and D fix the boundaries while ^„ and m„ govern the dynamics 
of the soliton interaction. In the case ctq = ^1 if c > and > the solution is regular. 

The masses of the solitons can be explicitly computed. For a generic choice of the 
parameters they do not depend on the spectral parameters A^,, p-n and on the initial 
positions (uom,Won). We get (ctq = -1) 

c I ^2 - f 



M. 



(-) 



M. 



(-) 



4 In [^1 + 4— 
41n ( 1 - 



M- 



(+) 
11 



(-) 
12 



4 In \\ 
4 In ^1 
41n II 



^1^2 ~ miTO2 p +4:cd I m2 p 
d I ^1 -Cm2 P 



£1^2 ~ mim2 p +4cc? | m2 
d I ^2 -Cmi P 



^1^2 



TOir7l2 



c \ £1 — Vmi 



-Acd I mi p 



£1^2 — miTO2 p +4cc? I mi p, 
4dcm2 + (^2^1 - m2mi)[(^i - Cm2)V - (mi - ^2)] 



£i^2 - mim2 p +4c(| m2 - £2!? p +d \ I2 P)] 
1 



M^i ' = 41n I 1 + 4dc— 



M 



(-) 
21 

(+) 
12 



^1^2 - mim2 P +4d(| m2 - liC P +c | £1 P)] 

|2 



m2 



^1^2 — mim2 
2 



= 4 In 1 + 4dc^ 



mi 



M^+^ = 41n ( 1 



I lit2 — mim2 
4dcmi + (4Z1 - m2mi)[(Z2 - Cmi)!? - (to2 - C^i) 



£i^2 - mim2 
1 



-4c(| mi -£il? P +d I li P)] 



[I 4^2 " mim2 P +4d(| mi ~ I2C P +c I ^2 P)] 



(3.159) 
(3.160) 
(3.161) 
(3.162) 

X 

(3.163) 
(3.164) 
(3.165) 

X 

(3.166) 
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Of special interest are the cases in which one or more masses are zero. For definiteness, 
we choose the masses M^nn {m ^ n) different from zero and we consider the case in 
which some M^nm are zero. 

For £i = Cto2 or £2 = 2?^2 one gets = or ' = and the solution 

describes the creation of a soliton. For £1 = Vrni or £2 — Cmi one gets M22^ = or 
M^^'' — 0, i.e. the annihilation of a soliton. For £1 — Cm2 and £2 — I?m2 it results 
that M22'' = M^^ •* = and the solution describes the creation of a pair of solitons. For 
£1 — Tfrni and £2 = Cmi one gets Mj^"* = M^^'' — , i.e. the annihilation of a pair of 
solitons. For £1 = Cm2 and £2 = Cmi one gets Af22^'' = m[^^ — 0, i.e. a soliton changes 
its mass and velocity. For £1 = Vmi, £2 = [\ V p /C)mi, m2 = {V/C)rfii one gets 
Mg^^ = m['^^ — and the solution describes the creation of a soliton. For £1 — Vfrii, 
£2 ^ {\ C p /P)m2, m2 = iV/C)mi one gets A'42^ = M^P = and_thc solution 
describes the annihilation of a soliton. For £1 = £2 — Cm, mi — m2 — m, Cm = Dm all 
the masses M^nm are zero and one gets a solution describing two interacting solitons. 

All these dynamical behaviours can be obtained by choosing the same boundaries. 
Boundaries do not give any information on the dynamics of solitons, but fix only their 
kinematics, i.e. their possible locations in the plane and their velocities. 

Of special interest is the case in which two spectral data are equal, say fii = iJ-2- In 
general, this solution describes the interaction of two solitons. The masses M^nmim = 
1,2) in this degenerate case depend also on the initial position of the solitons. With a 
special choice of the parameters one can get a solution describing the fission of a soliton 
(see figure 2 in |8|) and the fusion of two solitons. If in addition one chooses det{pri) = 
the solution describes a single soliton that by the interaction with a virtual soliton is 
forced to change velocity (see figure 1 in j^). Note that the two dynamical processes 
described in figure 1 and figure 2 in |^ are obtained by choosing the same boundaries 
and different matrices p and 77. 



3.4 Asymptotic bifurcation of multidimensional solitons 

The degeneracy of the solution when two discrete spectral data are chosen to be equal is 
worth of a deeper analysis. 

This phenomenon is well known in one dimension, but, while in one dimension by 
taking two eigenvalues equal in the A^-soliton solution we recover the (TV— l)-soliton 
solution, in 2+1 dimensions, as shown by the previous example, we get a new solution. 

To understand the underlaying mechanism it is convenient to consider the simplest 
case in which the effect takes place. Precisely, we choose N = 2, the matrices p and 77 
diagonal 

P diag(pi,p2), 77 = diag(77i,772) (3.167) 

and 

L) = ]1, C = ]1. (3.168) 

In the reduced case ctq = — 1 (we arc considering) the two matrices rj and p are related 
by the constraint 

PnPnn = -Vn^nn- (3.169) 

The set of these solutions describes a family of geometrical objects evolving in time in 
the {x, y) plane; each object E(^p f) of the family is parameterized by a point p in a P 
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space of 16 real parameters and by the time and therefore corresponds to a point in a 17- 
dimensional space S = {p,t}. We call generically stable or generic those objects of 
the family that depend in a difFcrentiable way on the 17 parameters in a neighborhood of a 
point So . According to the usual definition in catastrophe theory we call the complement 
of this open set {sq} the set of bifurcation points. 

The generic solution describes two solitons mutually interacting without changing 
shape and velocity. The only effect of the interaction is a shift in the position and in the 
overall phase. 

The only bifurcation points arc {p, — oo) and {p, +oo) with a special choice of the 
parameters p. Precisely, when any couple of discrete eigenvalues A„ , or /x„ , /x^ have 
the same real part, i.e. when the parameters belong to the hyperplanes Xnfn = XmfR or 
Mn» = MmJf i'^ the two-soliton solution is not stable at large times. For this special 
choice of the parameters one gets solitons that, as a result of their mutual interaction, 
exhibit a two dimensional shift and also a change of form. If we require, in addition, 
to the representative point p of the two-soliton solution to belong to the hyperplanes 
of lower dimensions Xn = Xm or /z„ = /z^ the two solitons, because of their mutual 
interaction, not only are shifted in the plane and change their form but also exchange 
mass. Surprisingly enough, in both cases the relevance of the bifurcation effect depends 
on the relative initial position of the two solitons. 

In order to describe a generic soliton with parameters A, /z and 7 = jTjp we introduce, 
in agreement with the notation for the one soliton solution, the two variables 



the phase 



= - Uo) - Xq{v - vo) + 2(A9AsR - msMsr)* 
(f) = nuiu - Uo) + Xu{v - Vo) + (A| - A| + /z| - 



(3.170) 
(3.171) 



and the functions 



a(±) = 7exp(±ei), ^ (1 +7)exp(6), b^'^ = 7exp(-6) 

£'(a, 6) = a^+^ + a^"^ + 6^+^ + b^~''- (3.172) 

To any of these variables and functions we add a label (n) when we are using the param- 
eters Xn, jJin and 7„ = jrjnPn of the nth soliton. 
Then the two-soliton solution can be written as 



(3.173) 



where 



N = 



^r/i-Aii exp(i^(i)) [^12 1^1^*12 P&(2)^ +t«i2Mi2^2iA2i&[2)^ + 



h^l2| V'12/''l2a(^)'' + |t*12p^2lA2ia|2)'' 



+ (1^2) 



(3.174) 



iAi2nMi2p&Sr>S2")^ 



+ |Al2|>1 



'■121 



'(1) «(2) 



\l) "(2) 



+ 



1^12|^|M12| 



"(1) "(2) ^ "(1) "(2) ^ "(1) "(2) 



+ 
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"(1) "(2) 



„(+)?,(+) +?,(+)„(+) 
" "(1) "(2) "t" '^(1) "(2) 



"(1) "(2) 



^(1) "(2) 



^(+)5(-) 
^(1) ^^(2) 



(1) «(2) 



-Re [??ir?2-^ii-A22^2i-!t*i2exp(i((^(i) - (?i(2)))] 



with 



(3.175) 

(3.176) 



A*12 — Ml ~ M2) Ai2 — Ai — A2. 

We consider the discrete values A„ fixed, with for instance A2Sft > Xm, and we study the 
bifurcations of the solution at large time with respect to the parameters There is 
a bifurcation at /xisr = iJ,2St and at ni = fi2- We need therefore to compute separately 
the asymptotic behaviour of the two soUton solution in the two cases. In particular we 
study the asymptotic behaviour g|^''(u, v) of the second soliton at < = ^poo. Thanks to 
the symmetry of the two-soliton solution the asymptotic behaviour of the other soliton 
is simply obtained by exchanging in the formulae the labels 1 and 2. 
For /xisR = /i2SR we get 

l+CXp(-2(/)[,'))S(2) 



9(2) 



q'12) = -2%A2aexp[i(;i(2))]- 



-2r?2A29 exp[i(0(2)-(/)[2)^)] - 



1 + T^^exp(~*0('))ii;, 



1+71 



(2) 



D^2)+J^E^2)D^2) 



(3.177) 



with 



(2) 



2) f(2) 
' S2 



exp 



M29 ^ ^ 



(2) 



C12 



C12 = exp[2,uiQ(uoi - 1*02) 



^(2) 



(2) A2) 
' ?2 



In 



M12 



so 



^12 
In 



\12 



12 



n, , /'f(2) ^(2) ^(2) ^(2)\ 

arg7ii2Ttti2 
argAi2Ai2. 



»(2) 



The mass of the second soliton at t = ±00 is 



(±) _ . i„ 1 + 72 



M(^2) =41n 



72 



(3.178) 



(3.179) 



Therefore the position and the phase of the soliton are shifted, its shape is changed but 
it docs not exchange mass with the other soliton. Note that its shape depends on the 
relative initial position of the two solitons. 
For fii = fj.2 we get 



9(2) 



„(+) _ 

9(2) - 



-2r72A2a exp[t((j!)(2) -0(2?)] 



(2) 



1+Ci 



^ 1+72 

' 72 



aJ+)+a[-;+6|J+(l+ci2)6j2l)] (C[ 
-2??2A2Qexp[i(/)(2)] 



2) ,(2)(-) .(2) ,(2)(-) 



?0 



I S2 SO 



7l 1+72 
! ^2 1+71 



"(2) "'""(2) " 



5(2) 



(1+C12 



1+71 



5(2) 



(3.180) 
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and 



M 




(2) 



4 In 



4 In 



(1+72)(1+C2l) 
Cl2+72(1 + C2l) 

(1 + 72)(1 + 71 + C2i7i) 



(3.181) 



72(1 + 7i) + C2i7i(l + 72) 



Therefore the soliton during the interaction shifts its position, changes shape and ex- 
changes mass with the other sohton, while the total mass of the two solitons is conserved. 
It is worth noting that the shape and the energy of each soliton depend on the relative 
initial position of the two solitons. 

Drawings describing the peculiar behaviours of solitons described in this section can 
be found in 

Finally let us remark that also the boundary a*^^-'(u,i) bifurcates. In fact, while in 
the generic case a'^-* (m, t) describes in the (m, t) plane two infinite waves crossing at one 
point, for /iijf = /h2sr it describes two parallel infinite waves and for /xi — /Lt2 only one 
infinite wave. 
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